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Abstract. Beta-delayed proton emission may occur at very low rates in the 
decays of the light nuclei ^^Be and *B. This paper explores the potential physical 
significance of such decays, estimates their rates and reports on first attempts to 
detect them: an experiment at ISOLDE/CERN gives a branching ratio for ^^Be 
of (2.5 ± 2.5) ■ 10~^ and an experiment at JYFL a 95% confidence upper limit of 
2.6 ■ 10-5 for **B. 
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1. Introduction 

The lightest chemical element in which beta-delayed proton (/3p) emission has been 
observed so far is carbon, in the decay of ^C. This paper explores the possibility of 
seeing the decay mode also in beryllium and boron, more specifically in the decays of 
-'^-'^Be and ^B. These nuclei are the lightest (not counting the deuteron) one-neutron 
and one-proton halo nuclei, respectively, and the possible (3p decays are intimately 
connected to this structure. There are two reasons for this connection, the first 
being the factorization of the wavefunction into a halo and a core part that suggests 
"independent decays" of the two parts [T] which naturally leads to final states with 
a continuum proton; we shall discuss this model in the next section. The second 
connection is specific for /3~ decays where one can derive the general expression [5] 

Q^pi^Z) = (m„ - mnjc' ~ Sn{^Z) = 782 keV - Sn{^Z) (1) 

from which it is seen that the /3p decay only occurs for nuclei with very low neutron 
separation energy. Similar relations exist for beta-delayed deuteron and triton emission 
that are known to take place close to the neutron dripline [H E] . So far /3p decays 
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have only been observed following or EC decays. General recent reviews of these 
processes can be found in [31 13] . 

The following section discusses the characteristics of the proposed decays in more 
detail including the estimated order of magnitude of the decay rate. The expected 
branching ratios are quite low and some general considerations are made on how they 
could be measured. Sections 3 and 4 give details on first measurements of the decays 
in ^^Be and '^B from which upper limits on the branching ratios can be given. The 
final section gives our conclusions and the prospects of identifying the two decays. 



2. Expected rates 

A crucial feature of halo states is their intrinsic clustering so that (a major part of) 
the wavefunction will factorize into a halo part and a core part O |6] . Formally, the 
beta-decay can then be thought of in terms of separate decays of the core (c) and halo 
(h) parts fll 

Ofi{\c)\h)) ^ {Op\c))\h) + \c){Op\h)) (2) 

but the actual final states will of course in general be superpositions of the two 
terms on the right side e.g. to ensure they have proper isospin. If the states on the 
right hand side are close to being eigenstates for the final system one would expect 
a large resemblance between the decays of the core nucleus and the halo nucleus. 
Experimental confirmation of this idea may be found in the decays of ^^Be and ^''Be, 
which are quite analogous [7]. 



2.1. ^^Be(l3p) 

The most recent mass tables [8] give — 11509.2 ±0.5 keV for the decay of 
-'^-'^Be, which leaves 281 keV as the energy window for beta-delayed proton emission. 
Different assumptions on the decay mechanism can give widely different rates [S] . We 
shall therefore present several simplified models that lead to rather straightforward 
expressions for the decay rate, and start by considering sequential decays as described 
in R-matrix theory that can handle broad intermediate states. However, the classical 
R-matrix expressions may not give a sufficiently accurate description of the decay 
mechanism. At the other extreme are models that consider decays to proceed directly 
to continuum states, as has been assumed for most calculations of the related /3d 
decays [2 [3]. 

Allowed beta decay of ^^Be will populate 1/2+ and 3/2+ levels in ^^B, but since 
these would decay to ^°Be-|-p with s-wave and d-wave protons, respectively, only the 
1/2+ states can be expected to give a sizeable contribution. Out of the known states 
[TUl [TT| in the excitation energy range 10.0-12.5 MeV the only one that may have 
spin and parity 1/2+ is a state at 11.444 MeV with width 103 keV that has been 
seen to decay by emission of a particles, but where the spin is so far undetermined. 
The R-matrix expressions for /3-decay have been derived for several cases, see [12] and 
references therein, but we shall here only consider a crude approximation based on 
the single-channel single-level case in which the energy spectrum is: 



_r 9l.,^ r/(27r) P{E-Et) 

wtot K [E - EoY +T^/AP[Eo- Et) 



where w denote the decay rates for the proton branch and the total beta decay, Ti the 
experimental halfiife, qa and K the standard weak decay constants [3l|4], the beta 
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decay phase space factor, Bqt the reduced Gamow-TeUer matrix elenient squared, Eg 
and r the level energy and width, and P{E ~ Et) the proton penetrability {Et being 
the proton threshold energy). Several energy dependent terms in the denominator 
are here approximated as constants. To take into account that the level can decay 
both by emission of a particles and protons one should multiply by Tp/Ttot to get the 
intensity of the proton spectrum. This ratio is unknown, but may in the best case 
(if the proton width is close to a single-particle width) reach 1/2. The GT strength 
Bgt can at most be 3 (the value for a pure neutron to proton decay). Inserting these 
maximum values the total branching ratio becomes 1.0 • 10~^, but a realistic value 
could be lower by an order of magnitude or more. The predicted spectrum is given 
in figure [1] for the case where the intensity (the product Bcr^p/Ttot) is a factor 50 
below the maximum, i.e. where the branching ratio is 2.0 • 10~^. 

As an intermediate step we consider next the suggested modification of the R- 
matrix formalism containing the opposite time ordering |13j . i.e. where the neutron is 
emitted before beta decaying into a proton. As argued by Barker |14j , who presents a 
more detailed derivation analogous to the way photon emission is treated, this may be 
a simple way to effectively include the decays to continuum states into the R-matrix 
framework. In the limit where only the "emission before decay" process is included 
the decay rate becomes 

— - T,/,-.UQ - E)BaT-^^-^^ , (4) 

where Sn is the neutron separation energy in ^^Be and the parameter 7^ quantifies 
the neutron emission probability. For a decaying level it is related to the level width 
through r = 2Pj^ and it must be less than the single particle width (the Wigner 
limit, here 5.8 MeV). The total branching ratio is proportional to 7^ and is at most 
2.5 • 10~*. The predicted energy distribution is also plotted in figure [TJ 

Finally, the simple model for decays directly to continuum states [15] gives 

^ = Tu4-f^iQ - ^)^oAE)^^^-^ , (5) 
wtot K R 271'' [hey 

where m is the reduced mass of the outgoing proton and R is the channel radius 
used to evaluate the penetrability. The Gamow- Teller strength is now explicitly 
energy dependent since the final state wavefunction is a continuum state. The 
expression above results if one uses a plane wave approximation for the final state when 
calculating Bgt{E). If one instead uses a Coulomb wave the effect of the transmission 
through the Coulomb potential is included in Bgt{E) and one must replace the ratio 
P/R by the wavenumber k. The total branching ratio is in this latter case 2.3 ■ 10~^, 
the corresponding energy distribution is given in figure [1] More realistic calculations 
will use distorted waves so that the final state strong interactions are also included. A 
recent more sophisticated calculation (16J in a two-body potential model gives a broad 
energy spectrum peaking at 0.1-0.2 MeV and a branching ratio of 3.0 • 10^*. (This 
specific calculation includes both Fermi and Gamow- Teller contributions to the decay 
rate.) 

To summarize, the different models predict somewhat different shapes of the 
energy distribution and give a branching ratio that typically is a few times 10~^, but 
could range up to 10~^. To see the decay mode an experiment therefore needs to be 
sensitive down to the 10~^ level, while considerably more intensity will be needed to 
separate the different models. 
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Figure 1. Energy spectra of beta-delayed protons from the decay of ^^Be 
predicted in three different models presented in the text. Sequential decay through 
the 11.444 MeV state in ^^B (dashed line), emission before decay (doted line), 
decay to continuum (full line). 

2.2. ^B(Pp) 

Turning now to ®B two facts make the decay directly to continuum states less likely 
here: the smaller spatial extent of the halo in this nucleus and the fact that there is 
an excited state in *Be within the /3-decay window (actually the EC-decay window) 
that the decay may pass through, namely the 1+ T = 1 state at 17.640 keV. It is 
situated 385 keV above the proton threshold and is know to decay mainly by proton 
emission [17) . If equation ([2]) is appropriate, we can estimate the matrix element of 
the transition to be the same as for the ground state decay of ^Be into the ground 
state of ^lA (the larger amplitude of the wavefunction in ^B increases the rate by 
a factor 2.1 [TH] but this is compensated by a factor 1/2 from the isospin Clebsch- 
Gordan coefficient squared). Both transitions are electron capture decays for which 
the phase space goes as the Q- value squared, so by scaling the halflife of ''Be we obtain 
an estimate of the branching ratio of the /3p transition of 2.3 • 10~*. 

Direct support for this estimate can be found from the three-cluster calculations 
in [m that predict a Bgt to the 17.64 MeV state of 1.366 and 1.997 for two different 
potentials employed. This is in the same range as the Bgt of 1.83 for the ^Be ground 
state transition. 
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Table 1. Q-values and separation energies in ^^Be beta-delayed particle decays. 



X 


P 


n 


a 


t 


S^("B) (keV)'^ 


11228.5 ±0.4 


11454.12 ±0.16 


8664.1 ±0.4 


11223.5 ±0.4 


(keV)" 


280.7 ±0.3 


55.1 ±0.5 


2845.2 ±0.2 


285.7 ±0.2 



"Mass values from ref. [8] 



2.3. Experimental considerations 

Direct detection of these /3p decays is challenging since one has to look for particles of 
low energy, with a low branching ratio, and where the decay branch has to be identified 
e.g. through particle identification (PID) techniques. As an alternative one may try 
to detect the presence of the daughter nucleus or, more ambitiously, its growing-in. In 
the cases in question here the final nucleus is stable or longlived so one cannot make 
use of its decay, but must resort to other methods. One possibility is accelerator mass 
spectroscopy (AMS) as explained in the next section. 

A low branching ratio by itself does not prohibit direct detection. As an example, 
branching ratios well below 10~^° have been detected for cluster radioactivity [20]. The 
challenge is rather to identify the decay branch uniquely. For proton kinetic energies 
below 1 MeV one could employ gas telescopes and use the energy loss AE versus full 
energy E to identify the particle. However, the energy loss for a proton peaks just 
below 100 keV energy so the AE-E curves for protons, deuterons and tritons will cross 
in this region, making separation e.g. between protons and tritons (relevant for ^^Be) 
impossible. Identification via detection of E and time-of-flight is an alternative, but 
will typically have much less efficiency. Separating protons from the background of 
energy loss signals from beta particles could also be a problem. 

If the proton can be recorded in coincidence with the recoiling nucleus this may 
be used to discriminate against other channels. This "ratio cut" method is useful to 
separate true events from response tails of events at higher energy, as done in the /3a 
decay of ^^N [3T] and needed here for the case of ^B. It can also be used to single out 
one decay channel among many others, as done for the ^"'^Li beta-decay to t+^Li [32] 
and needed here for ^^Be. The main challenge with this method is to obtain efficient 
detection of the recoil with sufficient accuracy in energy, in the examples given the 
energy was higher than for our /3p decays. 

3. The ^^Be experiment 

Four beta-delayed particle branches are energetically open for ^^Be, their Q-values are 
listed in table [TJ We may neglect the /3n branch since it would proceed to the ^°B 
3"*" ground state and require a d-wave neutron with very low energy. The /3t branch 
leads to ^Be and would therefore give a three-body final state, this could also be a 
quite interesting decay mode. The main background for the /3p branch will be /3a 
whose branching ratio was measured to be 3.1 ± 0.4% [23] (a recent determination 
gave 3.47 ±0.12% 
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3.1. The set-up and source production 

Due to the challenges involved in direct detection of the /3p decays we decided 
to explore the possibility of indirect detection. The /3p-daughter nucleus, ^°Be, is 
radioactive with a halfiife of 1.39 ■ 10^ years and is only present on earth in minute 
quantities due to production by cosmic rays. It is used as a tracer in earth sciences 
and procedures have been developed for detecting very small samples of ^'^Be via 
accelerator mass spectroscopy (AMS), see e.g. [25l[26]. Stated briefly, our experiment 
then consists in collecting a large number of -'^^Be atoms, measure their number by 
detecting 7-radiation from the ^^Be decay and later measure the number of produced 
^°Be atoms via AMS. 

The -'^^Be source was collected in September 2001 at the ISOLDE facility at 
CERN. The ^^Be atoms were produced by proton bombardement in a Ta target, 
ionized by laser ionization [27j . accelerated and transported to the measuring station 
where the ion beam passed two sets of collimators before being implanted into a Be 
foil (Goodfellow, LS226536 L C, thickness 0.01mm, purity 99.8+%). The average 
intensity of the collected beam was 3.6 • 10^ ions/s, the collection took place for close 
to 40 hours. The distance from the first collimator to the collection point was 177 
mm. The transmission was optimized with a stable ^^Na beam and a radioactive ^^Mg 
beam, the fraction of the activity that deposited on the collimator was at most 12%. 

3.2. Source strength 

Most of the decays of ^^Be proceed through states in ^^B that decay by gamma 
emission, the most prominent 7-ray being the 2124.47 keV line with a total branching 
ratio, including feeding from higher levels, of bj — 0.355 ± 0.018 [23]. The gamma 
detection was done with a Ge-detector that was placed about 40 cm from the collection 
point in the opposite direction of the collimator. It was further shielded by 3 mm Pb 
to reduce the total count rate. Still, the multi channel analyzer used for data taking 
had a deadtime of 14 ± 1% determined from the ratio of live time and real time of the 
data taking system. 

The Ge detector was energy and efficiency calibrated in situ with absolutely 
calibrated sources of ^°Co, ^"^^Cs and ^^^Th that bracket the energies from 239 keV to 
2615 keV. The total gamma spectrum from ^^Be is shown in figure [2] The spectrum 
is dominated by gamma lines from the decay of ^^Be with a small background mainly 
from "'"K and the ^^^Rn decay chain. The 2124 keV line is clearly separated from 
background and gives, after correction for detector efficiency and deadtime, a deduced 
amount of ^^Be atoms collected of (5.2 ± 0.3) • 10^^. As a cross-check of this number 
we shall look at two other 7-lines recorded from the ^^Be decay. First the 2895 keV 
line that again is nicely separated in our spectra; we deduce an intensity ratio with 
respect to the 2124 keV line of 0.241 ± 0.011% in agreement with the literature value 
of 0.227 ±0.008% [ig. Secondly, we consider the 478 keV line that follows /3a decays 
of -'^^Be to the first excited state in ^Li. The line is therefore recoil broadened and 
care must be taken when extracting its intensity since it furthermore is situated on 
top of a large background. Our deduced intensity relative to the 2124 keV transition 
is 0.75 ± 0.06%. 

The literature value for this last ratio can be derived from the decay scheme in 
[2H] to be 1.13% (it was not observed directly since gamma detection was done with 
a Nal detector). However, as we shall argue now, this value is most likely wrong. 
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Figure 2. Gamma energy spectra recorded from the decay of ^^Be, the number 
of counts per 0.4 keV is shown versus gamma-ray energy. The star marks the 478 
keV line. 



Alternative determinations of this ratio measured with Ge detectors can be extracted 
from experiments on the decay of ^^Li in which ^^Be is fed. From two such recent 
experiments accurate resuhs can be obtained (ahhough not pubUshed previously) : an 
experiment we did at ISOLDE [29] gave 0.67 ±0.11% and an experiment at TRIUMF 
[30j|3T] gave 0.69 ±0.07%. An unpublished recent f3a experiment [24^ has revised the 
branching ratio of alpha particles to the 478 keV excited state in ^Li to be 8.0 ± 0.4% 
of all alpha decays, a result based on a detailed analysis of the alpha spectrum but in 
disagreement with the literature value of 12.6 ± 1.2% [55]. Using the new value one 
derives a 478/2125 ratio of 0.78 ± 0.06%. Since the three new determinations all are 
consistent with our value for the ratio, we conclude that our spectra are internally 
consistent which adds confidence to our derived source intensity. 

3.3. Source purity 

The collection foil must of course be free of previous ^°Be activity, but we also need to 
worry about contamination coming along with the ^^Be beam. The most worrisome 
contaminant is ^^Li that can be ionized by surface ionization in the hot cavity where 
laser ionization of Be takes place. The primary production yield is lower by a factor 
at least 2000, but the main decay branch of ^"'^Li is beta-delayed one-neutron emission 
that leads to ^°Be and suppression of ^^Li is therefore needed. 

The mass difference between "Be and "Li, AM = Q("Li) = 20.551 MeV [8], 
is sufficiently big that the resolving power of ISOLDE'S high resolution separator 
HRS, which is up to M/AM = 5000, can easily separate the two nuclei that have 
M/AM — 500. A tail of ^^Li activity may nevertheless be present at the ^^Be mass 
position, but we have estimated its magnitude by looking for the corresponding tails 
of the nuclei ^Li and ^Li a factor M{^^Be)/M{}^Li) below their nominal setting. 
In both cases background prevented a direct identification — there are no 7-rays in 
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Table 2. Results of AMS measurements. 



Subfoil 


Mass (mg) 




i°Bc (10*^ atoms/mg) 


A (blank) 


3.04 


0.00149(12) 


2.67(22) 


B (background) 


3.65 


0.00152(11) 


2.73(20) 


C (catcher) 


4.22 


0.00169(13) 


3.03(23) 



either decay, so only the (non-unique) /3-detection is possible — but a lower limit 
of the suppression factor of 10'^ was extracted in both cases. A second source was 
collected with the separator deliberately set off-mass by two thirds of the ^^Li-^^Be 
mass difference (below the ^^Be position). The intensity of the ^^Be gamma lines was 
here reduced by a factor 500, consistent with the results at the lower masses. 

A further reduction can be obtained by blocking collection for the first 150 ms 
after the radioactivities are produced by the proton pulse. Most ^^Li, of halflife 8.5 
ms, will have decayed by then. However, at the end of our collection it was discovered 
that the onset of the blocking by mistake was delayed by a few ms allowing up to 
15% of the ^^Li to leak through. Thus up to 2 • 10^ ^°Be could have contaminated the 
sample through decay of collected ^^Li. 

A second possible contaminant could be the BeH molecule ^H^°Be that would 
be almost on top of the ^^Be mass since Af/AAf = 36600 for this case. It is quite 
unlikely that this molecule will be formed and ionized in sufficient quantities to be a 
problem, but the possibility should be kept in mind in future experiments. As will be 
seen we do not have indications for any of these contaminants in the present sample. 

3.4- AMS measurement 

The ""^"^Be accelerator mass spectrometry (AMS) measurements were performed at the 
Tandem Laboratory, Uppsala University, in Sweden |32] . A dedicated setup based on a 
NEC 5 MV tandem pelletron accelerator devoted to high precision and low background 
detection of ^'^Be was employed. The measurements were carried out according to 
the principle presented by Middletone and Klein [33l |34] where ion currents were 
normalized by -'^''O from the ^Be^'^O molecule. The absolute transmission in the AMS- 
system was determined with the NIST SRM 4325 standard (i°Be/^Be= 3.06 • 10"" 
[53] : NIST certificates the ratio to be Rsrm = 2.68 • 10~"). 

The Be-foil used in the investigation was divided into three 10 x 10 mm^ 
samples where one subfoil A was just chemically prepared and used as a blank value 
representative for the whole AMS procedure. Subfoil B and C were both mounted in 
two aluminium frames with a circular aperture (diameter 8 mm) . Both were installed 
in the ISOLDE facility. There was no implantation in subfoil B while subfoil C was 
used as the catcher for ^^Be. 

All three subfoils were chemically prepared by dissolution in HCl followed by 
precipitation to Be(0H)2 by adding NH3 which in turn was converted to BeO. Finally 
the BeO was mixed with Nb as a binder before the AMS-analysis in the accelerator. 

A summary of the results is given in table [2] No significant difference is observed 
for the blank and background foils which indicates no additional ^"^Be contamination 
in the handling in ISOLDE and during transportation. It is, however, noteworthy 
to compare the background values for the foil used with what is obtained from our 
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normal Be carrier which for 1 mg sample gives a R = 0.0005(1) corresponding to a 
factor two lower ^"^Be content. 

If we use subfoil B as a background to subfoil C the ^"Be signal can be estimated 
to be (1.28 ± 1.29) ■ 10^ atoms. With a total number of ^^Be nuclei in the sample of 
(5.2 ± 0.3) • 10^^ this leads to a deduced branching ratio of (2.5 ± 2.5) • 10~^ where 
the uncertainty is on the one sigma level. There is therefore no positive evidence for 
the ^^Be /3p decay, only an upper limit of the order of the most optimistic theoretical 
branching ratios. 

4. The ®B experiment 

The *B experiment was carried out at the IGISOL (Ion Guide Isotope Separator On- 
Line) separator at the Accelerator Laboratory of University of Jyvaskyla. The main 
aim of this experiment was to determine precisely the neutrino spectrum from the 
decay of ®B through analysis of the beta-delayed alpha spectrum. The result of this 
part of the experiment along with a detailed description of the experimental procedures 
is described in [36l[37]; the set-up consisted of four DSSSD detectors each backed by 
a thick Si detector in close geometry around the collection point. As a byproduct 
the high statistics of the experiment allowed to identify for the first time an electron 
capture branch in the decay of ®B, namely to the 2+ 16.922 MeV state. We are here 
concerned with the possible electron capture branch to the even higher-lying 17.640 
MeV state. 

Since the ^B experiment recorded a total of 16 million decays, the expected 
number of proton detections is 0.36 and we can only report an upper limit. The 
17.640 keV state will decay by emitting a 337 keV proton and a 48 keV recoiling ''Li 
ion that could not be detected. The main background from the decay will be a-a 
coincidences at low energy or positrons that happen to deposit around 340 keV in the 
detector. The latter background can be suppressed by requiring anti-coincidence with 
the backing detectors. The former background needs to be addressed indirectly since 
the set-up did not include particle identification. One can suppress the a particles by 
doing anti-coincidence with the opposite DSSSD, but to do this efficiently one needs 
to restrict the solid angle and only consider events close to the center of each DSSSD 
since beta-recoil can cause the a coincidences to deviate from a linear geometry. This 
reduces the effective solid angle by a factor about 9. A small background nevertheless 
remains, see figure[31 and the 95% confidence level upper limit of the number of counts 
in a 337 keV peak turns out to be 18. As seen in the figure there are indications for 
a peak at 349 keV; however, the uncertainty on the energy scale is only 7 keV since 
there were accurately known calibration lines from the decay of ^^Al close by |38) . 
The 95% confidence level upper limit on the number of counts in the 349 keV peak is 
45. Combing all factors gives an upper limit of 2.6 • 10~^ on the /3p branching ratio 
at 95% confidence level. To reach the range of predicted intensity the experiment 
needs to improve by a factor 1000. One may gain a factor of ten by including 
particle identification, which will reduce the background in the interesting energy 
range. Another factor of ten may be obtainable by increasing the effective solid angle, 
but one still would need a factor of ten higher overall production of ^B ions. 
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Figure 3. The low-energy spectrum from the decay of *B after background 
suppression. The red curve shows the two sigma upper limit from of a fit where 
the peak position was allowed to vary within uncertainties. 

5. Outlook 

The theoretical estimates in section 2 showed that beta-delayed proton emission in 
the decays of ^^Be and *B is likely to take place with a branching ratio that can be 
expected to be a few times 10~®. The energies of the emitted proton and the recoiling 
nucleus will be small. These two factors make detection of the processes challenging. 

We have carried out first experiments to look for these decay branches and have 
obtained upper limits on the branching ratios. These do not yet reach the level of the 
theoretical predictions, improvements are needed by factors of at least 10^ and 10^ 
for the two cases. The continuing improvements of radioactive beam facilities may 
give increases in the yields that will allow experiments to reach the needed level of 
sensitivity, but other possible improvements should also be considered. As explained 
above for changes to a more dedicated set-up should be considered. For the case 
of ^^Be the AMS measurement can be considerably improved if the amount of stable 
Be in the sample is reduced. This can be achieved by using a different catcher foil 
(e.g. Cu, Au or Nb) and after adding a small amount of Be carrier (< 0.1 mg with 
low ^'^Be content) do a refined chemical separation [35,. An enhanced sensitivity of 
two orders of magnitude can be expected from such improved sample handling. 

Once the decay modes have been established and the branching ratio securely 
determined, the next important step will be to determine the energy distribution of 
the emitted protons. This will give a more sensitive test of the theoretical calculations 
and may in particular show whether our hypothesis of a direct relation between the 
/3p decay mode and the halo structure in these two nuclei is correct. 
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